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Asstract. The paper presents the results of molecular analysis of aphids morphologically 
identical with Forda rotunda, a species hitherto not recorded in Polish fauna. The molecular 
analysis indicates that DNA sequences of “F. rotunda” (mitochondrial gene cox/ and nuclear 
efla) were identical with sequences of Mediterranean populations of Paracletus cimiciformis. 
The analysis thus suggests presence of a morphologically different form of P. cimiciformis in 
Poland. The applied molecular analysis proved its usefulness for solving taxonomic problems 
in a highly polymorphic group like aphids. 
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INTRODUCTION 


Paracletus cimiciformis von HEYDEN, 1837, a representative of tribe Fordini of 
family Eriosomatidae, is a species of European range. In the Mediterranean zone, on 
the territory of its primary host — Pistacia terebinthus — it is holocyclic and diecious 
(NieTo-NarriA et al. 2003). Bisexual generation reproduces and feeds on pistacia, 
while exules feed on different species of grasses, where they reproduce exclusively 
parthenogenetically. In Europe north of the Mediterranean, only such a permanently 
anholocyclic populations, feeding on grass roots, occur (SZELEGIEWICZ 1968). 

In these conditions these aphids are strictly myrmecophilous — they occur only in 
nests of Tetramorium caespitum (L.) and in Spain also in nests of the ant genus Tapinoma 
(NiETo-Narria et al. 2003). The morphs inhabiting ant nests are easily distinguished, 
with very characteristic 6-segmented antennae and white color of body, which is also 
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strongly flattened. Research suggests that this morph does not need to feed on plants 
as it can obtain food from ants on the way of trophallaxis (ZWOLFER 1958, HE 1980). 
All hitherto reported observations of this species in Poland concern the presence of this 
morph (SzELEGIEWwIcz 1968, PLacuta et al. 1996, Ostapacz & Haras 2009). 


1. The morphology of the studied form (‘A’) of the Paracletus cimiciformis 
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During the authors’ research on the anholocyclic populations of the tribe Fordini, 
we noticed that despite extensive field studies including a few dozens of sites with 
xerothermic vegetation, we found no specimens of P. cimiciformis accordant to the 
above description of its anholocyclic morph. 

Instead, the specimens with morphological features of the genus Forda were 
recorded, which should be determined as Forda rotunda Theobald, 1914, an aphid 
species never recorded in Poland. Those aphids were already reported in different parts 
of Europe, inhabiting nests of T. caespitum and feeding on the grass roots. Collected 
apterous viviparous females were round and shining green, with 5-segmented antennae 
(Ortiz-Rivas et al. 2009) (Fig. 1.). 

However, Roserti (1983, 1993) and later OrtTIz-RIvAs ET AL. (2009), basing on the 
molecular studies of Iberian representatives of the genus Paracletus, suggested that F 
rotunda may be a morph (so called form “A”) of P. cimiciformis. Because the factor 
triggering the development of the green morph is unknown, the molecular analysis 
was necessary to resolve the taxonomic identity of the studied aphids. For this pur- 
pose, two molecular markers were sequenced: mitochondrial cox/ and nuclear efla. 
These sequences were then compared to sequences of the taxon of interest deposited 
in GenBank by other authors. 

Mitochondrial genes have been widely used in the molecular systematics of in- 
sects (Simon et al. 1994). Due to their maternally inherited characteristics (Hoy, 1994), 
mitochondrial sequences are suitable for reconstructing phylogenetic relationships of 
parthenogenetic organisms like aphids (von DoHLEN & Moran 2000). Mitochondrial 
genes encoding cytochrome oxidase subunit 1 (cox/) and subunit 2 (cox2) have been 


Table 1. Collection data and accession numbers of used sequences. 


A : GenBank Accession No. 
. nt species 
k Voucher . N Collection Host g 
Species | Form Collection site associated 
number date plant ik aphid 
with aphids gla axl 
2 ` Poa Tetramorium 
= AA Dla Lubaczow 23.09.2010 . JQ241785 JQ241782 
5s Sp. caespitum 
= 
aS 
& E Poa Tetramorium 
a AA Dib Lubaczów 23.09.2010 . JQ241786 JQ241783 
2 Sp. caespitum 
3S 
v 
8 Poa Tetramorium 
Š AA D3 Sandomierz 22.09.2010 ; JQ241787 JQ241784 
Sp. caespitum 
Paracletus cimiciformis FM163597.1 AY227089.1 
Forda formicaria DQ499608.1 AY227086.1 
Species Forda marginata OSMANI EU701668.1 
ie P EU701669.1 
Forda riccobonii AY227076 1 
Pemphigus bursarius DQ779154.1 EU701834.1 
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used in phylogenetic studies of the three aphid subfamilies, Hormaphidinae (STERN 
1994; 1998; STERN et al. 1997), Lachninae (NoRMARK 1999; 2000) and Aphidinae (von 
Donen et al. 2006). These genes appear to be very useful for clarifying phylogenetic 
relationships also at the tribal and generic levels in aphids (NorMARK 2000; STERN 
1994). Among nuclear markers, also the nuclear gene encoding elongation factor |-alfa 
(efla) has often been used for species identification and phylogenetic reconstruction 
(Ortiz-Rivas et al. 2009). 

The aim of the study was to check the species identity of collected specimens of 
F. rotunda using molecular markers. 


MATERIAL AND METHODS 


The material comprised the three specimens of F rotunda (which were the supposed 


form “A” of P. cimiciformis), collected from two different sites located in southeastern 
Poland (Table 1). 


MOLECULAR RESEARCH 


Aphids were collected in 95% - 100% ethanol and stored at 4°C. We extracted 
total genomic DNA using the DNeasy Tissue Kit (Qiagen, Hilden, Germany). Prior 
to extraction, we removed the entire abdominal tagma to avoid possible contamination 
from gut content. 

Voucher specimens were collected in 80% ethanol and deposited in Zoology De- 
partment of the University of Silesia. Samples for extraction consisted of individuals 
from the same colony and were probably members of the same clone as the voucher 
specimens. 

Polymerase chain reaction (PCR) products were generated from a mitochondrial 
gene: cytochrome oxidase I (cox/), and nuclear gene: elongation factor | -alfa (ef/a). The 
primers used for the PCR amplifications were applied following published sequences 
(listed in Table 2), for cox/ using primers LCO1490 and HCO2198. The mitochondrial 
gene fragment was amplified as follows: initial denaturation at 94°C for 1 min, followed 


Table 2. The sequences of primers applied in the studied species 


. Gene Primers : a 
Region . Primer sequence (5 - 3’) Reference 
location name 


cox] |LCO1490 5'-GGTCAACAAATCATAAAGATATTGG-3' Folmer et al. 1994 
HCO2198 5'-TAAACTTCAGGGTGACCAAAAAATCA-3' 


Mitochondria 


Nuclear DNA fI Efs175 5'GGAAATGGGAAAAGGCTCCTTYAAGTAYGCYTGGG-3' Moran et al. (1999) 
efla 
` efrl 5'-GTGTGGCAATSCAANACNGGAGT-3’ Ortiz Rivas et al. (2009) 
7 
Nuclear DNA ae Efs1 nested 5'-TGGACAAAY TKAAGGCTGAACG-3' Ortiz Rivas et al. (2009) 


ted 
MIO) pet nested 5!-GTRTASCCRTTGGAAATTTGACC-3' 
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by 35 cycles of 94°C for 30s; annealing temperatures 47 for 45s; extension at 72°C for 
Imin; final extension at 72°C for 3 min. 

Two sets of nested primers were used for ef/a amplification. The first PCR was 
made using primers Efs175 and Efrl with the following conditions: initial denaturation 
at 94°C for 3 min, followed by 35 cycles of 94°C for 30s; annealing temperatures 52 
for 1 min; extension at 72°C for 2 min; final extension at 72°C for 7 min. Second, 
nested PCR was necessary, using primers Efs1, Efr3 on 0,51 taken from the first PCR. 
Conditions were identical except that the annealing temperature was 56°C. 

Products were sequenced at the Laboratory of DNA Sequencing and Oligonucle- 
otide Synthesis (IBB PAN). 

Amplified DNA was visualized using 1.5 % agarose gel electrophoresis with 
ethidium bromide staining. PCR products were purified using a QIAquick® PCR 
purification Kit (QIAGEN) and sequenced directly with an automated sequencer (Ge- 
nome Sequencer GS FLX Roche) at IBB PAN (Laboratory of DNA Sequencing and 
Oligonucleotide Synthesis). Complete sequences were deposited in GenBank under 
accession numbers given in Table 1, together with the accession numbers of sequences 
used for data comparison. 


D1B Paracletus cimiciformis 
e D3 Paracletus cimiciformis 
at DIA Paracletus cimiciformis 
gi|37955029|gb|AY227089.1| Paracletus cimiciformis 


nn) gi|188503307|gb|EU701669.1| Forda riccobonii 
gi|37954997|gb|AY227076.1| Forda riccobonii 

65 SS gil37955023|gb|AY227086.1| Forda formicaria 
93 gi|188503305|gb|EU701668.1| Forda marginata 


gi|188503637|gb|EU701834.1| Pemphigus bursarius 


2. Maximum parsimony (MP) tree from the analysis of the cox1 gene sequences from species in the genera 

Forda and Paracletus. Bootstrap values (as a percentage) are indicated when higher than 50. There were a 

total of425 positions in the final dataset, out of which 49 were parsimony informative. Phylogenetic analyses 
were conducted in MEGA4 


se ~ D1A Paracletus cimiciformis 


100 = DIB Paracletus cimiciformis 


D3 Paracletus cimiciformis 


73 ——————————— gi|208007561|emb|FM163597.1| Paracletus cimiciformis 
gi|98975523|gb|DQ499608.1| Forda formicaria 
99 ~~~  i|208007559|emb|FM163596.1| Forda marginata 
100 —— gi|98975525|gb|DQ499609.1| Forda marginata 


gil111013859|gb|DQ779154.1| Pemphigus bursarius 


3. Maximum parsimony (MP) tree from the analysis of the efla gene sequences from species in the genera 

Forda and Paracletus. Bootstrap values (as a percentage) are indicated when higher than 50. There were a 

total of 907 positions in the final dataset, out of which 64 were parsimony informative. Phylogenetic analyses 
were conducted in MEGA4 
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COMPUTER ANALYSIS OF DNA SEQUENCES 

The target sequences were 425 bp from the mitochondrial gene encoding cyto- 
chrome oxidase subunit (cox/), and 907 bp of the nuclear gene encoding elongation 
factor 1-alfa (ef7a) including both exons and introns. 

Chromatograms were analysed and assembled using the CHROMAS V2.3 
(TECHNELYsIUM Pty Ltp. 2004). Multiple alignments were done with Clustal X v1.81 
(THompson et al. 1997) with gap opening and gap extension penalties of 10.0 and 0.2, 
respectively, and subsequently manually revised. 


MOLECULAR PHYLOGENETIC ANALYSIS 

The phylogenetic trees were constructed using MEGA4: Molecular Evolutionary 
Genetics Analysis (MEGA) software version 4.0 (Tamura et al. 2007). 

Aligned sequences were tested using bootstrap method and to construct the tree the 
maximum parsimony (MP) method was applied. The MP tree was obtained using the 
close-neighbor-interchange algorithm (NEI & Kumar 2000) with search level 3 (FELSEN- 
STEIN 1985) in which the initial trees were obtained with the random addition of sequences 
(100 replicates). The codon positions included were 1 st+2nd+3rd+Noncoding. 


RESULTS 


The results of molecular analysis obtained with both applied markers indicated 
that sequences of studied specimens create a clade with those deposited in the data 
bank belonging to P. cimiciformis, which was supported with high bootstrap values. 
The results are presented on generated phylogenetic trees (Figs. 2 and 3). The clade 
consisting of a P. cimiciformis sequences is a sister clade to the one comprising the 
representatives of the genus Forda. The sequence of Pemphigus bursarius was applied 
as an outgroup as a representative of the sister subfamily Pemphiginae. 


DISCUSSION 


The sequences of analyzed specimens of the P cimiciformis form a clade in case 
of both applied markers, together with the sequences obtained from GenBank. This 
evidence seems sufficient to be sure that the form ‘A’ is a morph of the P. cimiciformis, 
and not a different species of the genus Forda. The Forda species constitute a separate 
clade, a sister group to the clade with Paracletus. Thus, the analysis of Central-European 
representatives of the tribe Fordini confirms the idea of synonymy of Forda rotunda 
with P. cimiciformis proposed by Ortiz-Rvas et al. (2009) on the basis of molecular 
research of Iberian populations. 

This is an interesting case of usefulness of molecular analysis for taxonomical pur- 
poses in aphids. The taxonomic problems of the species identity usually concern very 
closely related species, which are often morphologically very similar and thus, making 
it difficult to unambiguously determine the species (WANG & Qiao 2009, RAKAUSKAS 
et al. 2011). In the studied case there are two forms which may easily be determined 
as belonging to two different (although closely related) genera. In such case, when the 
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polymorphism ofa single species reaches the genus-level, only molecular analysis may 
correctly confirm the species identity. 

The results also clearly indicate that the bionomics of this species is not well rec- 
ognized yet. The factors determining the origin of both forms are still unknown and 
their ecological function is still a matter for further research. 
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